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FOREWORD 


The  model  study  reported  herein  was  requested  by  the  U.  S.  Army 
Engineer  Division,  North  Pacific,  in  a  letter  to  the  Office,  Chief  of 
Engineers ,  dated  20  April  19 66  and  was  subsequently  approved  in  a  letter 
to  the  North  Pacific  Division  dated  22  July  1966.  Authority  to  initiate 
the  investigation  was  granted  by  the  U.  S.  Army  Engineer  District, 
Seattle,  in  a  letter  to  the  Director,  U.  S.  Army  Engineer  Waterways  Ex¬ 
periment  Station  (WES) ,  dated  2k  January  1968. 

Design  and  construction  of  the  model  were  accomplished  during  the 
period  February  1968-June  1969;  hydraulic  and  salinity  verifications  were 
carried  out  during  the  period  June  1969-December  1969;  and  fixed-bed 
shoaling  verification  was  carried  out  during  the  period  August  1970- 
September  1970.  Hydraulic,  salinity,  and  dye  dispersion  base  tests  were 
performed  during  the  period  January  1970-July  1970.  After  completion  of 
air  fixed-bed  model  verifications  and  base  tests,  the  general  investiga¬ 
tion  program  was  initiated. 

This  report  describes  the  problems  that  necessitated  the  model  in¬ 
vestigation,  the  model  and  its  appurtenances,  and  the  fixed-bed  model 
verification  and  base  tests.  Subsequent  reports  will  describe  the  vari¬ 
ous  studies  conducted  in  the  model. 

The  study  was  conducted  in  the  Hydraulics  Division  of  the  WES 
under  the  direction  of  Mr.  E.  P.  Fortson,  Jr.,  Chief  of  the  Hydraulics 
Division,  and  Mr.  H.  B.  Simmons,  Chief  of  the  Estuaries  Branch.  The 
design,  construction,  verification,  and  base  tests  were  conducted  by 
Mr.  N.  J.  Brogdon,  Jr.,  Project  Engineer,  and  technicians  of  the  Estu¬ 
aries  Branch,  Tinder  the  supervision  of  Mr.  F.  A.  Herrmann,  Jr.,  Chief 
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of  the  Harbor  Entrance  Section.  This  report  was  prepared  by  Mr.  Brogdon 
with  the  assistance  of  Mr.  Herrmann. 

Personnel  of  the  Seattle  District  who  were  principally  responsible 
for  planning  the  course  of  study  included:  Messrs.  Dwain  Hogan*  John 
Norman,  Bruce  McCartney,  Bill  Peterson,  and  Harry  Disbrow.  These  per¬ 
sonnel  and  others  participated  in  conferences  on  the  model  study  and 
visited  the  model  from  time  to  time. 

Directors  of  WES  during  the  conduct  of  the  study  and  the  prepara¬ 
tion  and  publication  of  this  report  were  COL  Levi  A.  Brown,  CE,  and 
COL  Ernest  D.  Peixotto,  CE.  Technical  Directors  were  Messrs.  J.  B. 
Tiffany  and  F.  R.  Brown. 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  KJ«?<SUREMENT 


3»  -  Jsh  units  of  measurement 

. ic  units  as  Milovs: 

Multiply 

u^ed  in  this  report 

I 

By 

can  he  converted  to. 

1 

To  Obtain 

inches 

2.54 

centimeter  si 

feet 

0. 304*5 ' 

meters 

miles  (U.  S.  statute) 

1.609344 

kilometers  . 

square  feet 

0.092903 

Bquare  meters 

square  miles  /’T  S.  statute) 

2.58999 

■square  kilometers 

cubic  yards 

0.764555 

'  1 

cubic  meters 

feet  per  second 

0.3048 

mecers  per  second 

cubic  feet  per  second 

0.02831685 

cubic-  meters  per  second 

1 

1 

» 

4 

1 

! 

viX 


i' 


i 


i 


w 


I 


I 


I  ! 

I 

SUMMARY 


Grays  Harbor  model,  a  fixed-bed  model  with  provisions  for  future 
donversion  to  a  movable-bed  model,  was  constructed  to  scales  of  1:50C 
horizontally  and  1:100  vertically  and  reproduced  all  of  Grays  Harbor, 
the  Chehalis  River  ,t;o  the  head  of  tided  influence  (South  Montesano) ,  and 
a  portion  of  the  Pacific  Ocean  adjacent  to  the  harbor  entrance.  The 
primary  purposes  of  the  overall  model  study  will  be  to  determine  the  ef¬ 
fects  of  rehabilitation  of  both  the  north  and  south  jetties,  investigate 
the  stability  of  Point  Chehalis ,  determine  the  effects  of  enlarging  and 
realigning  the  navigation  channel,  and  locate  suitable  dredge  spoil  dis¬ 
posal  areas.  The  flushing  rate  and  dispersion  characteristics  of  waste 
materials  discharged  into  the  system  will  also  be  investigated. 

The  model  verification  tests  described  herein  indicated  that  the 
model  hydraulic  and  salinity  regimens  were  in  satisfactory  agreement 
with  those  of- the  prototype  for  comparable  conditions.  Iu  therefore  can 
be  assumed  that  the  model  will  provide  quantitative  answers  concerning 
the  effects  of  the  proposed  improvement  plans  on  the  hydraulic  and  sa¬ 
linity  regimens  of  the  estuary. 
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Fig.  1.  Vicinity  map 


GRAYS  HARBOR  ESTUARY ,  WASHINGTON 
VERIFICATION  AND  BASE  TESTS 

Hydraulic  Model  Investigation 


PART  I:  INTRODUCTION 
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The  Prototype 

1.  Grays  Harbor,  Wash.  (fig.  l),  located  about  45  miles*  north  of 
the  mouth  of  the  Columbia  River  and  93  miles  south  of  Cape  Flattery,  is 
a  large  bay  at  the  head  of  which  is  the  Chehalis  River.  The  roughly 
pear-shaped  harbor  diverges  from  the  Chehalis  River  at  Aberdeen ,  Wash . , 
15  miles  east  of  the  entrance,  to  a  maximum  width  of  13  miles,  including 
North  and  South  Bays  (plate  1).  The  water-surface  area  of  the  harbor 
varies  from  about  94  square  miles  at  mean  higher  high  water  (mhhv)  to 
about  35  square  miles  at  mean  lower  low  water  (mllw) .  The  harbor  is 
separated  from  the  Pacific  Ocean  by  two  narrow  sandy  peninsulas .  The 
northerly  peninsula  is  7  miles  long  and  terminates  at  Point  Brown,  Just 
north  of  the  harbor  entrance;  Point  Chehalis  comprises  the  nortn  end  of 
the  southerly  peninsula,  which  is  4  miles  long. 

2.  The  entrance  to  Grays  Harbor  (between  Point  Browr  and  Foint 
Chehalis)  is  about  2.5  miles  wide  and  is  protected  by  two  converging 
rubble-mound  Jetties  that  are  6500  ft  apart  at  their  outer  ends.  The 
existing  navigation  project  for  Grays  Harbor  and  Chehalis  River  was  au¬ 
thorized  by  the  Rivers  and  Harbors  Acts  of  3  June  1896,  2  March  1907, 

25  June  1910,  8  August  1917,  21  January  1927,  3  July  1930,  30  August 
1935 ,  2  March  1945,  30  June  1948,  and  3  September  1954.  The  project 
provides  for  a  600-ft-wide  by  30-ft-deep  entrance  channel  from  the 
Pacific  Ocean  across  the  bar  to  be  secured  by  dredging,  by  a  south  jetty 
13,700  ft  long,  and  by  a  north  jetty  16,000  ft  long.  The  top  elevation 


*  A  table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  is  presented  on  page  vii. 
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of  both  jetties,  as  originally  constructed,  was  +16  ft  mllw.  The  en¬ 
trance  channel  is  actually  self-maintaining  for  the  authorized  dimen¬ 
sions  and  thus  requires  no  maintenance  dredging.  Maximum  depths  in  the 
channel  are  about  -70  ft  mllw.  Based  on  the  recommendations  of  the 
U.  S.  Army  Engineer  Committee  on  Tidal  Hydraulics  and  the  results  of  a 
hydraulic  model  study  conducted  at  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  during  the  period  1950-1952,  it  was  concluded  that 
the  outer  6000  ft  of  the  south  jetty  should  be  allowed  to  deteriorate  to 
about  elevation  0.0  ft  mllw  to  alleviate  the  very  serious  erosion  prob¬ 
lem  at  Point  Chehalis  which  resulted  from  the  reconstruction  of  this 
jetty  during  the  period  1935-1939.  The  above-mentioned  Rivers  and 
Harbors  Acts  further  provide  for  maintenance  of  a  channel  30  ft  deep  and 
350  ft  wide  from  the  entrance  to  Cow  Point  (immediately  east  of  Rennie 
Island) ,  a  distance  of  14  miles ;  thence  30  ft  deep  and  200  ft  wide ,  for 
a  distance  of  4-1/8  miles;  a  turning  basin  30  ft  deep,  550  ft  wide,  and 
1000  ft  long  at  the  upstream  end  of  the  30-ft  channel  near  Cosmopolis ; 
thence  a  channel  16  ft  deep  and  150  ft  wide  for  10-3/4  miles  to  Monte- 
sano,  with  a  turning  basin  at  Montesano.  The  interior  channel  below 
Cow  Point  requires  annual  maintenance  dredging  of  about  1,600,000  c-u  yd. 

3.  Tides  occuring  in  Grays  Harbor  are  of  the  mixed  type  and  vary 
in  mean  diurnal  range  from  9.0  ft  at  Point  Chehalis  to  9.9  ft  at  Aberdeen. 
At  Point  Chehalis  the  elevation  of  mllw  is  -4.9  ft  mean  sea  level  (msl), 
while  that  of  mhhw  is  +4.1  ft  msl.  The  estuary  is  partially  mixed,  with 
a  significant  difference  between  surface  and  bottom  salinity  concentra¬ 
tions  but  no  distinct  saltwater  wedge.  Saltwater  intrusion  extends  up¬ 
stream  from  South  Aberdeen  (20  miles  from  the  entrance)  at  normal  fresh¬ 
water  flows  with  the  salinity  concentration  at  s’oa  18  (plate  l)  varying 
from  C.O  ppt  (total  salt)  at  lower  low  water  (llw)  to  about  10.0  to 
12.0  ppt  at  higher  high  water  (hhw).  The  horizontal  density  gradient 
from  the  entrance  to  the  upstream  limits  is  fairly  uniform.  Surface 
salinities  in  the  entrance  area  average  about  1.0  to  2.0  ppt  lower  than 
those  at  the  bottom,  while  near  the  upstream  limits  of  saltwater  intru¬ 
sion,  the  surface  salinities  average  about  3.0  to  5.0  ppt  lower  than 
those  at  the  bottom. 
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Purpose  of  Investigation 


In  general,  the  primary  purposes  of  the  overall  model  study 
are  to:  determine  the  need  for  and  develop  optimum  plans  for  rehabilita¬ 
tion  of  existing  jetties,  in  conjunction  with  erosion  studies  of  Point 
Chehalis  and  scour  adjacent  to  the  south  jetty;  determine  the  effects  of 
realignment  of  the  entrance  and  Sand  Island  shoal  navigation  channels; 
determine  the  effects  of  enlarging  the  navigation  channel;  determine  the 
location  of  suitable  dredge  spoil  disposal  areas;  investigate  ways  of 
reducing  maintenance  dredging  requirements  in  the  navigation  channel; 
determine  the  effects  of  enlarging  existing  small-boat  basins;  and  de¬ 
termine  the  flushing  rate  and  dispersion  characteristics  of  waste  mate¬ 
rials  that  are  discharged  into  the  estuary. 

5.  Following  the  rehabilitation  of  the  north  and  south  jetties  in 
19^*2,  the  entrance  channel  migrated  south  to  its  present  .^cation  adja¬ 
cent  to  the  south  jetty.  Scour  adjacent  to  this  jetty  has  caused  fear 
that  the  structure  may  eventually  be  undermined  unless  the  scour  can  be 
controlled  or  a  new  entrance  can  be  developed  in  the  northern  portion  of 
the  entrance.  The  alignments  of  the  entrance  and  the  bar  channels,  and 
their  close  proximity  to  the  south  jetty,  present  a  navigation  problem 
for  small  and  large  craft  alike  when  entering  and  leaving  the  harbor. 

This  problem  could  also  be  resolved  if  a  new  entrance  channel  were  de¬ 
veloped.  The  erosion  problem  at  Point  Chehalis  was  also  investigated  in 
conjunction  with  the  jetty  and  entrance  channel  studies. 

6.  Model  studies  will  be  conducted  to  develop  plans  to  reduce  the 
existing  annual  maintenan. a  dredging  requirements  in  the  interior  channel 
of  the  harbor.  These  studies  will  be  conducted  in  conjunction  with  pro¬ 
posed  changes  in  channel  alignments  and  enlarged  channel  dimensions. 
Studies  will  be  made  for  both  existing  and  proposed  channel  changes  to 
determine  best  locations  of  spoil  disposal  areas  and  to  determine  if 
turning  basins  along  the  channel  might  also  serve  as  sediment  traps. 

This  report  contains  all  important  data  related  to  the  hydraulic,  salin¬ 
ity,  and  fixed-bed  shoaling  (entrance  area)  verification  of  the  model. 
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PART  II:  THE  MODEL 


Description 

7.  The  Grays  Harbor  model  reproduced  approximately  230  square 
miles  of  the  prototype  area,  including  the  Chehalis  River  to  South  Mon- 
tesano,  the  Pacific  coast  from  about  7-5  miles  north  of  the  north  Jetty 
to  about  7.5  miles  south  of  the  south  jetty,  and  offshore  areas  veil  be¬ 
yond  the  -60  ft  contours ;  Elk,  Johns,  Humptulips ,  Hoquiam,  and  Wishkah 
Rivers ;  North  and  South  Bays ;  pod  the  system  of  sloughs  and  creeks  that 
affects  tidal  action  throughout  the  model  area.  The  limits  of  the  area 
reproduced  are  shown  in  plate  1.  A  general  view  of  the  model  is  shown 
in  fig.  2. 

8.  The  model  was  constructed  to  linear  scale  ratios,  model  to  pro¬ 
totype,  of  1:500  horizontally  and  1:100  vertically.  From  these  basic 
ratios  the  following  scale  relations  were  computed  by  the  Freudian  rela¬ 
tions:  slope  5:1,  velocity  1:10,  time  1:50,  discharge  1:500,000,  and 
volume  1:25,000,000.  The  salinity  and  dye  concentration  ratios  for  the 
study  were  1:1.  One  prototype  cycle  (diurnal  tide)  of  2h  hr  and  50  min 
was  reproduced  in  the  model  in  2 9  min  and  1*8.5  sec.  Horizontal  grid  co¬ 
ordinates  are  based  on  Polyconio  Projection  1927  North  American  Datum, 
and  vertical  control  was  based  on  USC&GS  msl  1929  adjustment.  The  model 
was  approximately  320  ft  long  and  180  ft  wide  at  its  widest  point,  and 
covered  an  area  of  about  30,000  sq  ft.  It  was  completely  enclosed  to 
protect  it  and  its  appurtenances  from  the  weather  and  to  permit  unintei- 
rup+ed  operation. 

9.  The  model  vas  initially  constructed  as  a  fixed-bed  model;  how¬ 
ever,  provisions  were  made  to  convert  the  entrance  area  to  a  movable -bed 
model  at  a  later  date  when  such  studies  were  deemed  necessary.  Limits 
of  the  movable-bed  section  are  shown  in  plate  1.  The  model  was  molded 
to  conform  to  the  prototype  hydrographic  conditions  that  existed  in  1967. 
The  navigation  channel  was  molded  in  removable  blocks  so  that  desired 
alterations  could  readily  be  made. 

10.  The  permanent  model  roughness  employed  consisted  of 
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General  view  of  model,  looking  upstream 


3A-in.-wi.de  metal  strips  placed  in  depths  greater  than  6  ft  below  msl 
and  cut  off  below  the  low-water  elevation.  The  use  of  these  metal 
strips  as  roughness  was  necessary  because  proper  adjustment  of  velocity 
and  distribution  of  currents,  both  horizontally  and  vertically,  in  any 
given  cross  section  could  not  be  obtained  by  the  use  of  ordinary  bound¬ 
ary  roughness  alone  in  the  deep  areas  of  the  model.  The  areas  above 
-6  ft  msl  (shoal  areas  and  tide  flats)  were  roughened  by  raking  the  model 
surface  during  construction  to  provide  the  desired  degree  of  roughness. 

Appurtenances 

11.  The  model  was  equipped  with  the  necessary  appurtenances  to 
reproduce  and  measure  all  pertinent  phenomena  such  as  tidal  elevations, 
saltwater  intrusion,  current  velocities,  freshwater  inflow,  waves,  dis¬ 
persion  characteristics,  and  shoaling  distribution.  Apparatus  used  in 
connection  with  the  reproduction  and  measurement  of  these  phenomena  in¬ 
clude  a  tide  generator  and  recorder,  tide  gages,  salinity  meters,  salin¬ 
ity  samplers,  chemical  titration  equipment,  current  velocity  meters, 
freshwater  measuring  weirs,  wave  generators,  dye  injection  and  measuring 
equipment,  and  shoaling  injection  and  recovery  apparatus.  This  equip¬ 
ment  is  described  in  detail  in  subsequent  paragraphs. 

Tide  generator  and  recorder 

12.  The  reproduction  of  tidal  action  in  the  model  was  accom¬ 
plished  by  means  of  a  tide  generator  (figs.  3  and  4),  located  in  the 
model  ocean,  that  maintained  a  differential  between  a  pumped  inflow  of 
salt  water  to  the  model  and  a  gravity  return  flow  to  the  supply  sump  as 
required  to  reproduce  all  characteristics  of  the  prototype  tides  at  the 
control  station  (Westport).  The  tide  generator  was  equipped  with  a  con¬ 
tinuous  tide  recorder  so  that  the  accuracy  of  the  model  tide  reproduc¬ 
tion  could  be  checked  visually  at  any  time. 

Tide  gages 

13.  Permanently  mounted  point  gages  were  installed  in  the  model 
at  the  locations  of  the  seven  recording  tide  gages  used  for  collection 
of  field  tide  data  (plate  l).  Portable  point  gages  were  used  to  measure 
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OPERATION  OF  TIDE  GENERATOR 

THE  WATER  SURFACE  IN  THE  MODEL  (A'  IS  APPRECIABLY  HIGHER  THAN  i  THE  SUMP  (B)  A  PUMP  (C)  LOCATED 
IN  THE  SUMP  DISCHARGES  A  CONSTANT  AMOUNT  OF  WATER  INTO  ONE  SI -C  01  A  OIVIDEO  HEAOBAY  (01.  DISCHARGE 
FROM  THE  MODEL  TO  THE  SUMP  IS  CONTROLLED  BY  AN  AUTOMATIC  VALVE  ‘.E  j  OPERATION  OF  THE  AUTOMATIC  VALVE 
IS  OICTATEO  BY  THE  TIOE  CONTROL  (FI  LOCATED  IN  THE  MODEL.  WHEN  THE  TIDE  IS  RISINS  IN  THE  MODEL.  THE 
AUTOMATIC  VALVE  IS  ALMOST  CLOSED  30  THAT  MOST  OF  THE  WATER  REMAINS  IN  THE  MODEL.  WHEN  THE  TIDE  IS  FALL¬ 
ING  IN  THE  MODEL,  THE  VALVE  IS  ALMOST  OPEN  SO  THAT  ALL  THE  WATER  FROM  THE  PUMP,  PLUS  THE  DESIRED  AMOUNT 
OF  WATER  FROM  THE  MODEL  RETURNS  TO  THE  SUMP  THE  TIOE  CONTROL  MAINTAINS  THE  PROPER  VALVE  OPENING  AT  ALL 
TIMES  AS  REQUIRED  TO  REPRODUCE  ANY  OESIREO  TIOE  IN  THE  MODEL 


Fig.  3.  Operation  of  tide  generator 


Fig.  h.  Tide  recorder 


tidal  elevation  at  other  points  as  required.  The  model  gages  were  grad¬ 
uated  in  0.001  ft  (0.1  ft  prototype). 

Salinity  and  dye  samplers 

lh.  Salinity  and  dye  water  samples  were  drawn  into  collection 
vials  by  negative  pressure  from  a  vacuum  pump  connected  to  a  central  man¬ 
ifold,  which  in  turn  was  connected  to  tubes  running  to  each  sampling  lo¬ 
cation.  This  device  enabled  simultaneous  sampling  at  all  desired  depths 
at  all  sampling  stations  throughout  the  model.  Details  of  the  multidepth 
.‘■ampler  are  shown  in  fig.  5 • 


Fig.  5.  Multidepth  sampler 


Chemical  titration  equipment 

15.  This  method  of  determining  salinity  concentration  was  used 
primarily  to  determine  the  salinity  concentration  at  the  saltwater  source 
(sump).  The  equipment  consisted  of  a  graduated  burette  for  measuring  the 
volume  of  silver  nitrate,  pipettes  for  measuring  the  volume  of  samples 
used,  sample  jars  in  which  to  perform  the  titration,  a  supply  of  silver 
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nitrate,  and  a  quantity  of  potassium  chromate  for  use  as  an  end-point 
indicator  in  the  titration  process.  The  method  consisted  of  adding  a 
known  concentration  of  silver  nitrate  solution  to  a  known  volume  of  the 
model  salinity  sample;  the  amount  of  silver  nitrate  required  to  precip¬ 
itate  the  salt  contained  in  the  sample  was  then  converted  to  salinity  in 
parts  per  thousand  of  NaCl. 

Salinity  meters 

16.  All  salinity  concentrations  for  samples  taken  from  the  model 
were  determined  by  use  of  conductivity  cells  especially  built  and  cali¬ 
brated  for  this  purpose.  The  salinity  meter  assembly  is  shown  in  fig.  6 


Fig.  6.  Salinity  meter  assembly 
Current  velocity  meters 

IT.  Current  velocity  measurements  were  obtained  with  miniature 
Price-type  current  meters  (fig.  7).  The  five  meter  cups,  constructed  of 
a  light  plastic  or  metal  material,  were  approximately  O.OU  ft  (U  ft  pro¬ 
totype)  in  diameter  and  were  mounted  on  a  horizontal  wheel  0.09  ft  in 
diameter;  the  center  of  the  cups  was  0.05  ft  (5  ft  prototype)  from  the 


Fig.  7*  Miniature  Price-type  current  meter 

bottom  of  the  frame.  The  meters  were  calibrated  frequently  to  ensure 
accurate  operation  and  were  capable  of  measuring  actual  velocities  as 
low  as  0.03  fps  (0.3  fps  prototype). 

Freshwater  inflow  measuring  devices 

18.  All  rivers  with  significant  mean  freshwater  inflows  were 
equipped  with  a  constant  head  tank  and  Van  Leer  weir  for  precise  mea¬ 
surements  of  their  respective  flows.  The  inflows  of  streams  with  minor 
freshwater  inflows  were  combined  with  those  of  nearby  streams  of  signif¬ 
icant  inflow,  and  the  combined  inflow  wss  introduced  at  a  central  point. 
Individual  inflow  points  were  located  on  the  Chehalis,  Hoquiam,  and 
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Wishkah  Rivers.  Chenois  Creek  and  Grass  Creek  inflows  were  combined  with 
that  of  Humptulips  River.  Johns  River,  Newskah  Creek,  Charley  Creek, 
Stafford  Creek,  Indian  Creek,  Campbell  Creek,  Chapin  Creek,  O'Leary 
Creek,  and  Andrews  Creek  inflows  were  combined  with  that  of  Elk  River. 
Skimming  weir 

19.  The  mixed  salt  water  and  fresh  water  that  accumulated  in  the 
model  ocean  had  to  be  removed  in  order  to  maintain  a  constant  volume  and 
a  constant  source  salinity.  This  was  accomplished  by  means  of  skimming 
weirs  that  removed  a  quantity  of  mixed  water  from  the  surface  layer 
equal  to  the  freshwater  inflow  to  the  estuary.  Precise  measurement  of 
the  combined  discharge  from  the  skimming  weirs  was  made  by  use  of  a 

Vein  Leer  weir. 

Wave  generators 

20.  The  model  ocean  was  equipped  with  two  20-ft-long  wave  gener¬ 
ators  to  produce  the  effects  of  ocean  waves  on  the  transportation  and 
deposition  of  sediments.  The  wave  generators  were  of  the  plunger  type 
and  could  be  adjusted  to  produce  the  desired  wave  height  and  period. 

Dye  injection  and  measuring  equipment 

21.  Model  tests  were  made  to  determine  the  flushing  rate  and  dis¬ 
persion  characteristics  of  the  model.  A  given  weight  of  powdered  fluo¬ 
rescent  dye  was  thoroughly  mixed  with  a  known  volume  of  water  and  then 
stored  in  a  glass-sided  tank.  The  tanks  were  equipped  with  a  system  of 
valves  and  tubes  to  control  the  desired  discharge  at  the  injection  lo¬ 
cations.  Injection  was  made  at  two  locations  simultaneously  by  utiliz¬ 
ing  two  different  types  of  dyes.  Water  samples  were  collected  at  lo¬ 
cations  throughout  the  model  with  a  sampling  device  identical  with  the 
multidepth  salinity  sampler  described  in  paragraph  1^.  Concentrations 
of  the  samples  were  measured  by  means  of  a  Turner  Model  III  fluorometer 
(fig.  8). 

Shoaling  injection 
and  recovery  apparatus 

22.  Shoaling  was  reproduced  in  the  model  by  injecting  granulated 
(l/8  in.  by  1/8  in.)  nylon  particles  or  polystyrene  plastics.  The  spe¬ 
cific  weight  of  the  model  sediment  to  be  used  for  a  particular  problem 
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Fig.  3.  Turner  fluotfometer  i 

was  determined  by  trial  and  error,  and  varied  from  about  1.03  to  1.20. 

1 

The  injection  procedure  for  individual  problem  areas  varied;  however, 
the  normal  methods  were  to  place  the  material  by  hand  or  allow  a  slurry 
of  the  material  to  flow  into  the  model  through  a  perforated  trough.  At 
the  end  of  each  model  test,  the  shoaling  material,  deposited  within  the 
limits  of  the  navigation  channel  or  prescribed  problem  area  was  ire- 
covered  and  measured  volumetrically. 
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:  .  PART  III:  VERIFICATION  OF  THE  MODEL 

23.  :  Verification  of  the  Grays  Harbor  fixed-bed  model  was  accom¬ 
plished  in  three  phases:  (a)  hydraulic  verification,  which  ensured 
that  tidal-  elevations  and  times,  and  current  velocities  and  directions 
were  in  proper  agreement  with  the  prototype;  (b)  salinity  verification, 
which  ensured  that  salinity  phenomena  in  the  model  correspond  to  those 

l 

,  of  the  prototype  for  similar  conditions  of  tide,  ocean  salinity,  and 

I 

freshwater  inflow;  and  (c)  fixed-bed  shoaling  verification,  which  en¬ 
sured  acceptable  reproduction  of  prototype  shoaling  distribution. 
Fixed-bed  shoaling  verification  for  the  interior  30-  by  350-ft  naviga- 
tipn  channel  and  movable-bed  shoaling  verification  are  beyond  the  scope 
of  this  report;  however,  the . fixed-bed  shoaling  verification  for  the 
area  between  the  jetties  will  be  discussed.  The  need  for  movable-bed 
studies  had  not  teen  determined  at  the  time  of  preparation  and  publi¬ 
cation  of  this  report;  therefore,  the  movable-bed  verification  had  not 

'  s  I  I 

been  initiated.  i 

.  l 

i  2k.  The  accurate  reproduction  of  hydraulic  and  salinity  phe¬ 
nomena  in  an  estuary  model  is  an  important  phas$  in  the  preparation  of 
the  model  for  its  ultimate  use  in  evaluating  the  effects  of  proposed 
improvement  works.  Every  effort  was  made  to  obtain  a  comprehensive 

verification  of  till  pertinent5  phenomena.  Prototype  data  used  for  the 

.i  ! 

model  hydraulic  and  salinity  verification  were  collected  in  the  field 
by  U.  S.  Army  Engineer  District,  Seattle,  personnel  during  the  period 
k  October  to  17  October  1967.  These  data  were  furnished  to  Waterways 
Experiment  Station  engineers  in  the  form  of  curves  and  tables.  The 
entrance  shbaling  verification  (fixed  bed)  wis  based  on  soundings 
taken  in  September  1965  and ’November  1969. 

25.  This  report  contains  all  important  data  related  to  the  hy- 

1 

araulic,  salinity,  and  fixed-bed  shoaling  (entrance  area)  verification 
of  the  model,  -in  order  to  facilitate  reference  to  these  data  in  other 
phases  of  the  pverall  investigation.  , 
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Hydraulic  Verification 


Prototype  data 

2 6.  In  1967  the  Seattle  District  undertook  an  extensive  prototype 
metering  program  of  the  Grays  Harbor  estuary  in  order  to  obtain  data 
with  which  to  adjust  and  verify  the  Grays  Harbor  estuary  model.  The 
metering  program  was  carried  out  in  two  consecutive  6-day  periods  that 
were  selected  as  being  generally  representative  of  normal  spring  and 
neap  tide  conditions  on  the  basis  of  predicted  tides.  Eighteen  stations 
on  eight  ranges  (plate  1)  were  established  at  which  to  collect  current 
and  salinity  data.  Because  of  insufficient  personnel  and  equipment  and 
prohibitive  cost,  it  was  impossible  to  obtain  simultaneous  observations 
at  all  stations. 

Tidal  adjustment 

27.  The  objective  of  the  model  tidal  adjustment  was  to  obtain  an 
accurate  reproduction  of  prototype  tidal  elevations  and  tidal  phases 
throughout  the  model.  Prototype  tidal  data  from  seven  recordirg  tide 
gages  (plate  l)  were  availah  '  to  verify  the  accuracy  of  the  model  tidal 
adjustment.  These  gages  recorded  essentially  continuously  throughout 
the  prototype  data  collection  period  discussed  previously. 

28.  During  each  of  the  6-day  metering  periods  there  were  signifi¬ 
cant  variations  of  tidal  ranges  and  other  tidal  characteristics.  In 
order  to  avoid  the  time-consuming  and  expensive  procedures  of  adjusting 
the  model  to  reproduce  all  of  the  tides  observed  during  the  metering 
program,  it  was  decided  to  select  one  24.84-hr  (diurnal)  tide  for  each 
of  the  6-day  periods  that  approximated  the  tidal  condition  for  that  par¬ 
ticular  period  and  to  complete  the  adjustment  of  tides,  currents,  and 
salinities  throughout  the  model  for  only  the  two  tides  thus  selected. 

The  two  tides  chosen  were  6-7  October  1967  (spring  tide)  and  13-14  Octo¬ 
ber  1967  (neap  tide). 

29.  The  procedure  followed  was  to  adjust  the  tide  generator  in 
such  a  manner  that  the  tides  generated  in  the  model  ocean  would  cause 
an  accurate  reproduction  of  prototype  tides  at  Westport  (control)  tide 
gage,  then  to  adjust  the  model  roughness  until  prototype  tidal 
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elevations  and  times  were  reproduced  to  scale  throughout  the  model. 

30.  Comparisons  of  model  and  prototype  tidal  data  for  the  two 
tides  reproduced  in  the  model  are  presented  in  plates  2-9.  Plates  2-k 
show  tidal  elevations  for  the  spring  tide  condition  at  the  Westport,  Bay 
City,  Goose  Island,  Crossover  Channel,  Aberdeen  Port  Dock,  South  Aber¬ 
deen,  rcid  Montesano  tide  gages.  Low-  and  high-water  lunar  intervals; 
high-,  low-,  and  mean-water  levels;  and  range  of  tide  profiles  are  pre¬ 
sented  in  plate  5  for  the  spring  tide  condition.  The  maximum  discrep¬ 
ancy  in  tidal  range  was  in  the  order  of  0.5  ft  prototype  (0.005  ft 
model).  This  maximum  discrepancy  occurred  at  the  gage  located  at  Monte¬ 
sano.  This  gage  also  reflected  the  greatest  discrepancy  for  high-  and 
low-water  lunar  intervals  and  high-  and  low-water  levels.  This  gage  was 
located  very  near  the  upstream  model  limits  (plate  1)  and  the  freshwater 
discharge  point  for  the  Chehalis  River.  Although  tidal  effects  extend  a 
considerable  distance  upstream  from  this  gage  location  in  the  prototype, 
there  was  no  supplemental  provision  for  the  passage  of  tidal  flow  at  the 
upstream  model  limit.  Since  the  tidal  flow  in  this  area  is  quite  smell, 
and  since  this  area  is  a  considerable  distance  upstream  from  any  poten¬ 
tial  problem  areeis  that  might  have  been  subject  no  model  investigations, 
it  was  believed  that  provision  for  such  tided  flow  would  have  involved 
an  unnecessary  expense.  The  discrepancies  between  model  and  prototype 
tidal  data  at  the  Montesano  gage  are  a  direct  result  of  the  close  prox¬ 
imity  of  this  gage  to  the  upstream  model  limit.  Comparisons  of  model 
and  prototype  tidal  data  for  the  neap  tide  condition  are  presented  in 
plates  6-9.  The  maximum  discrepancy  in  tidal  range  was  about  1.2  ft  and 
was  observed  at  the  Montesano  gage  location. 

Adjustment  of  currents 

31.  The  objective  of  the  model  current  adjustment  was  to  obtain 
an  accurate  reproduction  of  prototype  current  velocities  and  distribu¬ 
tions  throughout  the  model.  Prototype  current  velocity  data  were  avail¬ 
able  at  18  stations,  the  locations  of  which  are  shown  in  plate  1.  Pro¬ 
totype  readings  were  made  at  the  surface ,  quarter  depth ,  middepth , 
three-quarter  depth,  and  2  ft  above  the  bottom  for  a  period  of  at  least 
25  hr  at  each  station.  As  discussed  above,  no  single  simultaneous 
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survey  of  the  entire  estuary  vas  made;  pud  the  model  adjustment  was 
made  for  only  one  tide  from  each  of  the  6-day  metering  periods.  The 
neap  and  spring  tides  selected  for  reproduction  in  the  model  corre¬ 
sponded  to  the  days  on  which  current  measurements  were  obtained  on 
range  3  (sta  7,  8,  9*  10,  11,  and  12).  Since  the  tidal  characteris¬ 
tics  varied  throughout  each  6-day  period  it  is  obvious  that  current 
velocities  observed  on  any  particular  day  are  not  necessarily  repre¬ 
sentative  of  those  which  would  have  been  observed  at  the  same  station 
on  any  other  day  within  the  6-day  observation  period.  Thus,  it  was 
necessary  to  adjust  at  least  the  maximum  prototype  velocity  observa¬ 
tions  in  order  that  they  would  more  nearly  represent  prototype  condi¬ 
tions  for  the  tides  reproduced  in  the  model.  It  was  found  that  there 
was  a  reasonably  good  linear  relation  between  tidal  range  and  maximum 
velocity.  The  adjusted  maximum  prototype  velocities  are  shown  in  the 
velocity  verification  plates  (10-73)  for  those  stations  at  which  the 
prototype  observations  were  made  on  a  different  tide  than  that  re¬ 
produced  in  the  model.  Also  shown  for  each  current  velocity  station 
are  the  tide  which  was  reproduced  in  the  model  and  that  which  actu¬ 
ally  occurred  during  prototype  data  collection  corresponding  to  the 
prototype  current  velocity  data  presented. 

32.  The  procedure  followed  for  adjustment  of  current  veloc¬ 
ities  was  to  reproduce  each  of  the  two  tidal  and  discharge  condi¬ 
tions  in  turn  and  adjust  the  model  roughness  until  the  current  veloc¬ 
ities  at  each  metering  station  were  correctly  reproduced  in  the  model. 
The  freshwater  discharges  used  during  the  current  velocity  verifica¬ 
tion  were  obtained  by  averaging  the  inflows  that  occurred  during  each 
of  che  6-day  metering  periods.  Comparisons  of  model  and  prototype 
current  velocities  for  all  stations  are  presented  in  plates  10-73. 
Measurements  obtained  at  half-hour  intervals  were  plotted  for  both 
model  and  prototype,  And  smooth  curves  were  drawn  through  the  points. 
No  attempt  will  be  made  to  discuss  each  comparison  of  prototype  and 
model  measurements,  but  the  agreement  obtained  throughout  the  model 
is  considered  to  be  very  satis factory. 


Salinity  Verification 


33.  The  objective  of  the  model  salinity  adjustment  was  to  obtain 
an  accurate  reproduction  of  vertical  and  lateral  distribution  of  proto¬ 
type  salinities  throughout  the  model.  Reproduction  of  prototype  salin¬ 
ity  phenomena  in  the  model  required  the  maintenance  of  the  proper  salin¬ 
ity  in  the  ocean  water-supply  system  and  the  establishment  of  the  proper 
mixing  environment.  Because  prototype  salinity  data  obtained  at  sta¬ 
tions  near  the  mouth  of  the  estuary  were  obviously  in  error,  it  was  de¬ 
cided  to  ’use  a  model  ocean  salinity  that  was  representative  of  salini¬ 
ties  known  to  *iccur  near  the  entrance  to  most  west  cosat  estuaries  (33.0 
ppt).  The  prototype  salinity  data  at  sta  1,  6,  and  10  for  the  neap  tide 
condition  and  those  at  sta  1  for  the  spring  tide  conditions  are  omitted 
from  this  report  because  they  were  obviously  erroneous.  Prototype  ob¬ 
servations  at  these  locations  were  all  made  with  the  same  meter.  Salin¬ 
ity  data  at  other  times  and  stations  obtained  with  the  same  meter  were 
silso  obviously  in  error,  but  they  were  not  used  in  the  model  verifica¬ 
tion.  Salinity  observations  were  made  at  hourly  (prototype)  intervals  in 
both  the  model  and  prototype.  These  data  were  plotted  and  smooth  curves 
were  drawn  through  the  points  and  are  compered  with  corresponding  proto¬ 
type  curves  in  plates  74-105. 

34.  The  agreement  demonstrated  between  model  and  prototype  is 
considered  excellent.  It  is  pointed  out  that  no  additional  adjustment 
of  the  model  was  necessary  to  obtain  the  agreement  shown  in  plates 
74-105.  This  substantiates  the  model  adjustment  of  tides  and  currents 
and  indicates  that  bhe  upland  fresh  water  w as  being  properly  mixed  with 
salt  water  from  the  ocean  supply. 

Entrance  Area  Fixed-Bed  Shoaling  Verification 

35.  The  entrance  area  fixed-bed  verification  involved  the  repro¬ 
duction  of  the  prototype  shoaling  pattern  throughout  the  entrance  area. 
The  basic  objective  of  the  model  shoaling  verification  was  to  identify 
a  synthetic  sediment  that  would  move  and  deposit  raider  the  influence  of 
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model  forces  in  the  same  manner  that  the  natural  sediments  move  and 
deposit  under  the  influence  of  natural  forces.  In  the  process  of  iden¬ 
tifying  a  suitable  model  sediment,  a  great  number  of  variables  were  in¬ 
volved  and  each  had  to  be  resolved  by  tried  and  error.  The  most  signif¬ 
icant  variables  included:  (a)  shape,  size,  gradation,  and  specific 
weight  of  the  artificial  sediment;  (b)  method,  location,  duration,  and 
quantity  of  artificial  sediment  injection;  (c)  height  and  direction  of 
waves;  (d)  magnitude  of  tide  range;  (e)  length  of  model  operation;  and 
(f)  readjustment  of  model  roughness.  The  final  adjustment  was  accom¬ 
plished  using  a  granulated  nylon  material  with  a  specific  weight  of 
1.13  to  1.15  and  a  grain  size  of  about  1/8  in.  by  1/8  in.  (cylinders). 

A  spring  tide  was  reproduced  in  the  model  with  a  mean  freshwater  dis¬ 
charge,  and  waves  were  generated  alternately  from  the  west  and  south¬ 
west.  Because  of  the  distortion  of  the  model  scales,  it  was  not 
possible  to  reproduce  waves  in  the  model  that  represented  prototype 
conditions.  At  the  start  of  the  shoaling  test,  100,000  cc  of  the  above 

material  was  injected  into  the  model  over  a  period  of  five  cycles  begin¬ 

ning  at  hour  22  of  cycle  1.  During  cycles  1  and  2,  5000  cc  was  placed 
in  the  model  at  hour  22  along  a  line  represented  by  the  east  boundari es 
of  shoaling  sections  18,  19,  C-7,  and  20  (plate  2);  lO.GGO  cc  was  placed 

in  the  model  along  this  same  line  at  hour  9-5;  and  5000  cc  was  placed  in 

the  model  along  a  line  running  north  and  south  through  the  middle  of 
section  19  at  hour  14.  During  cycles  3,  k ,  and  5,  5000  cc  was  placed  in 
the  model  at  hour  22  along  a  line  represented  by  the  east  boundaries  of 
shoaling  sections  18  and  19;  10,000  cc  was  placed  in  the  model  along  this 
same  line  at  hour  9*5,  and  5000  cc  was  placed  in  the  model  at  hour  14 
along  a  north-south  line  represented  by  the  center  line  of  section  19. 

It  was  necessary  to  operate  two  wave  machines  in  the  model  ocean  to 
simulate  the  effects  of  wave  energy  on  the  resuspension  and  movement  of 
the  model  sediment.  Wave  machine  1  was  oriented  so  that  waves  were  gen¬ 
erated  from  the  west;  wave  machine  2  was  oriented  so  that  waves  were 
generated  from  the  southwest.  The  wave  heights  and  periods  did  not  rep¬ 
resent  any  particular  prototype  waves ,  but  they  were  adjusted  to  simu¬ 
late  the  degree  of  agitation  required  so  the  model  shoal  materials  could 


be  moved  and  deposited  by  the  tided,  currents.  Wave  machine  1  vaa  oper¬ 
ated  between  hours  2  and  17;  wave  machine  2  was  operated  between  hours 
20  and  2.  This  wave  machire  operating  procedure  was  repeated  during 
each  of  the  five  injection  cycles.  At  the  end  of  the  test,  the  material 
was  recovered  from  individual  sections  (plate  106)  and  measured 
volumet  ri c  ally . 

3 6.  Scour  and  fill  patterns  that  developed  in  the  prototype  dur¬ 
ing  the  period  1965-1969  are  shown  in  plate  107.  Only  the  fill  patterns 
indicated  in  this  plate  were  used  as  a  basis  for  the  entrance  shoaling 
verification  since  scour  cannot  be  reproduced  in  a  fixed-bed  model.  The 
relative  shoaling  pattern  achieved  in  the  model  is  shown  in  plate  108. 
The  actual  volumes  of  material  retrieved  from  individual  sections  are 
shown  in  table  1  along  with  the  respective  percentages  of  the  total 
amount  injected  and  recovered. 


Limitations  of  the  Accuracy  of  Model  Measurements 


37*  Measurements  of  tidal  elevations  in  the  model  were  made  with 
point  gages  graduated  to  0.001  ft,  or  0.1  ft  prototype.  The  limitations 
of  the  current  velocity  meters  used  in  the  model  should  be  considered  in 
making  close  comparisons  between  model  and  prototype  velocity  data.  The 
center  line  of  the  meter  cup  was  about  0.05  ft  above  the  bottom  of  the 
frame;  therefore,  bottom  velocity  measurements  in  the  model  were  actu¬ 
ally  obtained  at  a  point  5.0  ft  (prototype)  above  the  bottom,  instead  of 
about  2.0  ft  as  in  the  prototype  metering  program.  The  model  veloci¬ 
ties  were  determined  by  counting  the  number  of  revolutions  in  a  10-sec 
interval  (which  represented  a  period  of  about  8  min  in  the  prototype), 
as  compared  with  about  a  1-min  observation  in  the  prototype.  The  hori¬ 
zontal  spread  of  the  entire  meter  cup  wheel  was  about  0.11  ft  in  the 
model,  representing  about  55  ft  in  the  prototype,  as  compared  with  less 
than  1.0  ft  for  the  prototype  meter.  Thus,  the  distortion  of  area 
(model  to  prototype)  results  in  comparison  of  prototype  point  velocities 
with  model  mean  velocities  for  a  much  larger  area.  The  same  is  true  for 
the  vertical  area,  since  the  height  of  the  meter  cup  was  about  0.0U  ft 
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(4.0  ft  prototype)  as  compared  with  only  a  few  Inches  for  the  prototype 
meter. 

38.  All  model  salinity  measurements  presented  in  this  report  were 
made  with  a  salinity  meter  (conductivity  type)  and  are  c<~  sidered  to  he 
accurate  within  0.5  ppt  in  the  higher  ranges  and  0.2  ppt  in  the  lower 
ranges.  The  model  samples  were  collected  at  the  bottom,  middepth,  and 
surface  elevations.  The  elevations  of  the  bottom  and  middepth  samplers 
were  fixed  in  the  model  and  were  not  allowed  to  vary  with  the  tide  as 
was  the  surface  sampler.  Simultaneous  water  samples  were  drawn  into 
vials  from  the  three  elevations  by  means  of  a  vacuum  system,  whereas 
the  prototype  salinities  were  measured  in  place  at  successive  depths. 
Similar  to  the  model  velocity  data,  the  model  salinity  data  also  repre¬ 
sent  an  average  over  a  much  larger  prototype  area,  Bince  the  vacuum  sam¬ 
pling  system  used  in  the  model  drew  the  sample  from  a  radius  of  about 
1/2  to  1  in.  (80  ft  in  the  prototype).  The  accuracy  with  which  the 
model  could  be  expected  to  duplicate  salinities  from  cycle  to  cycle  for 
identical  conditions  appears  to  be  about  +3  percent. 

Discussion  of  Results  of  Verification  Tests 


39.  Agreement  between  model  and  prototype  phenomena ,  as  evidenced 
by  the  results  of  hydraulic  and  salinity  verification  data,  appears  to 
be  excellent.  The  model  was  considered  to  be  sufficiently  similar  to 
its  prototype  to  be  confidently  utilized  in  quantitative  studies  of  the 
effects  of  proposed  improvement  plans  on  hydraulic  phenomena  and  salin¬ 
ity  intrusion. 

40.  In  fixed-bed  shoaling  tests  it  is  not  possible  to  reproduce 
bed  scour.  Should  any  significant  scour  areas  develop  in  the  prototype 
as  the  result  of  some  new  construction,  they  would  create  a  source  of 
new  sediments  that  could  cause  significantly  increased  shoaling  in  other 
areas  that  would  not  be  predicted  from  the  model  results.  The  accuracy 
of  the  nodel  shoaling  tests  is  considered  to  be  about  +10  percent,  since 
that  is  about  the  limit  of  accuracy  of  repeating  identical  tests. 


•  > 
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PART  IV:  BASE  TEST  RESULTS 


41.  In  order  to  evaluate  the  effects  of  the  proposed  Improvement 
plans,  it  was  first  necessary  to  establish  hydraulic,  salinity,  shoal¬ 
ing,  surface  current  pattern,  and  dye  dispersion  "base  tests”  that  de¬ 
picted,  respectively,  the  characteristics  of  these  conditions  throughout 
the  model  for  existing  conditions.  Thus,  a  test  in  which  no  improvement 
plan  was  installed  in  the  model  is  referred  to  as  a  base  test,  since  its 
results  constitute  a  basis  of  comparison  for  determining  the  effects  of 
improvement  plans.  The  model  base  tests  were  conducted  with  model  con¬ 
ditions  of  a  mean  tide,  a  source  salinity  concentration  of  33.0  ppt,  and 
a  total  freshwater  discharge  of  11 ,1+00  cfs .  The  discharges  at  individual 
inflow  points  were:  Chehalis  River,  6300  cfs;  Hoquiam  River,  970  cfs; 
Wislikah  River,  420  cfs;  Humptulips  River,  2750  cfs;  and  Elk  River,  960 
cfs.  Prior  to  data  collection  of  any  type,  the  model  was  operated  for 

a  sufficient  period  of  time  to  achieve  a  condition  of  dynamic  salinity 
stability.  The  procedure  for  starting  operation  of  a  tidal  model  uti¬ 
lizing  seat  water  is  developed  during  the  verification  phase  of  the 
study  and  thereafter  is  followed  exactly.  The  initial  salinity  concen¬ 
trations  are  artificial;  and  the  model  must  be  operated  for  a  period  of 
time  to  allow  salinities  to  reach  stable  values  with  respect  to  time, 
depth,  and  location.  It  was  found  for  the  Grays  Harbor  model  that  the 
best  procedure  was  to  flood  the  entire  model  with  ocean  water  to  the  el¬ 
evation  of  mhhw,  then  to  turn  on  the  freshwater  inflow  weirs  and  allow 
the  model  to  operate  until  salinity  stability  was  achieved.  It  was  nec¬ 
essary  to  operate  the  model  for  about  five  tidal  cycles  (about  2.5  hr) 
before  relatively  stable  conditions  existed,  after  which  data  collection 
could  be  initiated.  To  ensure  a  higher  degree  of  salinity  stability 
when  obtaining  salinity  measurements ,  salinity  samples  were  not  taken 
until  the  model  had  operated  for  at  least  10  cycles  (about  4  hr). 

42.  The  locations  of  all  stations  monitored  during  base  test  con¬ 
ditions  are  shown  in  plate  109.  Tide  gage  locations  are  identical  with 
those  in  the  verification,  except  that  two  gages  were  added  in  the  en¬ 
trance  area  (Ocean  and  Westport  "A”)  to  more  accurately  define  tidal 

•  ►*< 
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changes  in  this  area  resulting  from  any  tests  investigated.  The  verifi¬ 
cation  of  current  velocities  and  salinities  was  based  on  the  results  of 
data  obtained  at  18  locations  throughout  the  estuary  (see  plate  l).  For 
obvious  reasons ,  none  of  these  stations  could  be  located  along  the  navi¬ 
gation  channel  center  line.  The  coverage  was  not  sufficient,  partic¬ 
ularly  in  the  entrance  area,  to  give  a  complete  composite  picture  of  the 
effects  of  future  improvement  plans  to  be  tested;  therefore,  the  station 
location  and  numbering  scheme  used  during  the  verification  phase  of  the 
model  study  was  abandoned  during  the  testing  phase  for  a  sampling  scheme 
that  would  better  define  the  effects  of  proposed  plans.  These  sampling 
stations  were  numbered  19-58  in  order  to  avoid  confusion  with  the  veri¬ 
fication  stations  (1-18).  Several  stations  were  added  after  the  above 
scheme  was  adopted  and  are  postfixed  with  a  letter  to  the  number  of  the 
nearest  station. 

Tides 

43.  Base  test  tidal  height  data  were  obtained  at  nine  locations 
throughout  the  model.  The  data,  collected  at  half-hour  intervals,  were 
plotted  and  smooth  curves  were  drawn  through  the  points.  The  locations 
of  tide  stations  are  shown  in  plate  109  and  data  are  presented  in 
plates  110-112. 

Current  Velocities 

44.  Base  test  current  velocity  data  were  obtained  at  44  locations 
throughout  the  model.  Bata  were  collected  at  the  surface,  middepth,  and 
bottom  elevations  at  half-hour  intervals  at  all  locations  where  the 
depth  of  water  was  sufficient.  At  the  few  stations  that  were  located  in 
shallow  water,  middepth  measurements  were  not  made.  The  half-hour  mea¬ 
surements  were  plotted,  and  smooth  curves  were  drawn  through  the  points. 
The  locations  of  current  velocity  stations  are  shown  in  plate  109,  and 
the  current  velocity  data  are  presented  in  plates  113-156. 


Salinities 


1*5.  Base  test  salinity  data  were  obtained  at  hi  locations 
throughout  the  model.  Water  samples  were  collected  at  the  surface,  mid¬ 
depth,  and  bottom  elevations  at  hourly  intervals.  The  salinity  concen¬ 
trations  were  determined  with  a  salinity  meter  and  were  later  plotted 
and  smooth  curves  were  drawn  through  the  points.  The  locations  of  sa¬ 
linity  stations  are  shown  in  plate  109,  and  the  salinity  data  are  pre- 

|  sented  in  plates  157-197. 
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'  Entrance  Shoaling 

} 

t 

1*6.  The  final  results  for  the  verification  of  shoaling  in  the 
j  entrance  area,  discussed  in  paragraphs  35  and  36  of  this  report ,  also 

|  served  as  the  shoaling  base  test ,  since  there  were  no  plans  installed  in 

|  the  model  during  the  verification  tests.  The  location  of  the  shoal 

!  sections,  prototype  data  used  as  the  basis  for  verification,  and  final 

model  shoaling  verification  results  are  shown  in  plates  106,  107,  and 
108,  respectively.  The  volumes  of  shoaling  material  retrieved  from  in¬ 
dividual  sections  are  presented  in  table  1.  The  shoaling  test  technique 
developed  for  this  model  provided  only  qualitative  information  concern¬ 
ing  entrance  area  shoaling.  It  was  possible  to  identify  areas  of  heavy 
shoaling  and  to  compare  the  merits  of  various  plans  on  the  basis  of  the 
relative  quantities  of  shoaling  in  the  entrance  area;  however,  it  was 
not  possible  to  relate  quantities  of  model  shoaling  to  prototype 
quantities. 

Bye  Dispersion 

1*7.  Model  tests  were  made  to  determine  the  flushing  rate  and  dis¬ 
persion  characteristics  of  the  estuary  for  two  very  generalized  sources 
(one  upstream  and  one  at  the  entrance)  of  materials  entering  the  estu- 
.  ary.  Test  procedures  consisted  of  introducing  one  type  of  conservative 

j  fluorescent  dye  (Pontacyl  Brilliant  Pink)  at  a  point  near  Cosmopolis  and 

J 
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another  type  (Uranine)  near  the  outer  end  of  the  south  jetty.  Four 
liters  of  Pontacyl  Brilliant  Pink  dye  with  an  initial  concentration  of 
10,000,000  pph  and  the  density  of  fresh  water  of  1.0  was  injected  into 
the  model  near  Cosmopolis  at  the  middepth  elevation  over  a  complete 
tidal  cycle  beginning  at  hour  0  of  cycle  1.  Four  liters  of  Uranine  dye 
with  an  initial  concentration  of  10,000,000  pph  and  the  density  of  sea 
water  vas  injected  into  the  model  near  the  outer  end  of  the  south  Jetty 
at  the  bottom  of  the  channel  from  hour  20  (llw)  of  cycle  1  until  hour  2 
(lhw)  of  cycle  2.  Each  dye  mixture  was  contained  in  a  glass  standpipe 
and  was  introduced  at  a  constant  rate  throughout  the  injection  period 
through  a  1/4-in.  copper  tube  oriented  at  the  proper  elevation.  Prior 
to  initiating  dye  injection,  the  model  was  operated  for  four  tidal 
cycles  to  establish  salinity  stability  throughout  lie  model. 

1+8 .  The  locations  of  the  two  dye  release  points  are  shown  in 
plate  109 •  During  the  dye  dispersion  tests,  water  samples  were  with¬ 
drawn  from  the  model  for  subsequent  analysis  for  dye  concentration.  The 
samples  were  taken  at  4l  locations  with  the  multidepth  sampler  described 
in  paragraph  14.  Samples  were  obtained  at  surface,  middepth,  and  bottom 
elevations  at  the  times  of  hhw  slack  current  (hour  13)  and  llw  slack 
current  (hour  20)  over  a  period  of  16  tide'  cycles.  Concentrations  of 
the  samples  were  measured  by  means  of  a  Turner  Model  III  fluorometer. 

1+9.  Plots  of  dye  concentration  as  a  function  of  time  at  typical 
stations  for  the  Cosmopolis  dye  release  are  shown  in  plates  198-202. 
Smooth  curves  have  been  drawn  through  the  data  points.  Plots  of  dye 
concentration  profiles  along  the  navigation  channel  after  5  and  12  tidal 
cycles  for  the  Cosmopolis  dye  release  are  shown  in  plates  203  and  204, 
respectively.  Similar  plots  of  concentration  histories  and  profiles  for 
the  south  jetty  dye  release  are  presented  in  plates  205-211.  Table  2 
presents  the  dye  concentrations  for  high-  and  low-water  slacks  at  sur¬ 
face,  middepth,  and  bottom  elevations  at  all  stations  for  both  the  Cos¬ 
mopolis  and  south  jetty  dye  releases. 
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Current  Pattern  Mosaics 


50.  Surface  current  pattern  mosaics  were  made  of  the  entrance 
area  for  use  in  evaluating  proposed  channel  realignments,  effectiveness 
of  dikes,  groins,  etc.  The  mosaic  also  provides  a  means  for  current 
velocity  analysis  in  areas  too  shallow  for  measurenents  with  the  veloc¬ 
ity  meter.  The  mosaics  included  in  this  report  cover  only  the  entrance 
area  and  a  portion  of  the  inner  harhor;  however,  it  was  anticipated  that 
as  the  model  study  progressed  this  partial  mosaic  would  be  expanded  to 
include  the  entire  estuary. 

51.  The  mosaics  were  prepared  from  time-exposure  photographs  of 
confetti  floating  on  the  water  surface.  A  bright  light  was  flashed  im¬ 
mediately  before  the  camera  lens  was  closed,  resulting  in  a  bright  spot 
at  approximately  the  end  of  each  confetti  streak  which  indicates  the 
direction  of  flow.  Current  velocities  can  be  determined  from  the  photo¬ 
graphs  by  measuring  the  lengths  of  the  confetti  streaks  and  comparing 
the  lengths  with  the  velocity  scale  presented  in  each  mosaic.  Photo¬ 
graphs  were  taken  at  hourly  (prototype)  intervals  throughout  a  complete 
tidal  cycle  (24.84  hr).  Surface  current  pattern  mosaics  for  base  test 
conditions  are  presented  in  photos  1-25. 
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Table  1 

Entrance  Shoaling 
Base  Test  Conditions 


Section 

No. 

Volume 

Retrieved 

cc* 

of  Total 
Injeetpd 

Section 

No 

Volume 

Retrieved 

1 

(  Percent 

of  Total 

1 

10 

0.01 

0.01 

31 

1,070 

Retrieved 

1.28  ! 

Injected 

1.07 

C. 

1,125 

1.36 

1.12 

32 

15 

1  0.02 

0.02 

3 

220 

0.26 

0.22 

33 

1  0 

0.00  . 

oioo 

4 

1,230 

1.47 

1.23 

34 

'  2,175' 

2.62 

2.18 

5 

5,680 

6.83 

5.68 

:  35 

Q 

0.00 

OiOO 

6 

10 

0.01 

0.01 

36  ' 

’  *  0 

0.00 

,  0.00 

7 

20 

0.02 

0.02 

.  37 

'  1,810 

2.18 

i.81 

8 

790 

0.95 

0.79 

38 

3,890 

4.68 

1  3.89 

9 

3,750 

4.50 

3.75 

39 

5  ■ 

0.00 

0.00 

10 

1,930 

2.32 

1.93 

40 

5 

6.00 

o-.oo 

11 

3,595 

4.31 

3.59 

4l 

0 

.  o.oo! 

6.00 

12 

15 

0.02 

42 

0 

■  0.00 

0.00 

13 

3,960 

4.76 

3.96 

43 

6,570 

7.90 

6.57i 

14 

3,535 

4.26 

3.54 

44 

0  . 

c.oo 

0.00 

15 

5,985 

7.19 

5.98 

1  .  45 

0 

0.00 

0.00 

16 

3,440 

4.13 

3,44 

46 

:  •  0  1 

:  0.00 

0.00 

17 

1,465 

1.76 

1.46 

.47 

835 

1.00 

0.84 

18 

495 

0.59 

O.49  • 

48 

2,710 

3.27 

2.71 

19 

2,440 

2.94 

2.44 

49 

880 

1.06 

0.88 

20 

3,070 

3.69 

3.07 

C-l 

15 

J 

0.02 

0i02 

21 

25 

0.03 

0.03 

■  C-2 

0 

0.00  1 

0.00 

22 

25 

0.03 

O.03 

■  C-3 

.  0 

0.00 

p.oo 

23 

1,210 

1.46 

1.21 

C-4 

15 

0.02  : 

0.02  , 

24 

965 

1.16 

O.96  1 

.  C-5 

3,795 

•.4.56 

3.79 

25 

'  340 

0.40 

0.34 

C-6  : 

885 

1.05 

6.38 

26 

10 

0.01 

0.01 

C-7 

1 

1,110 

1.3.4 

•  1.11 

27 

2,435 

2.94 

2.44  , 

C -8 

3,035 

3.65 

1 

,  3.04 

28 

2,190 

2.64 

2.19 

C-9 

2,245 

2.69 

2.24 

29 

2.16 

1.80 

C-10 

340 

0.42  . 

0.1?4 

30 

15 

0.02 

0,02 

c-11 

15 

6.01 

0.01 

Total 

83,200 

109.00 

83.20  1 

Note:  Total  volume  injected  was  100,000  cc, 
*  Average  of  two  identical  runs. 
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VERIFICATION  OF 
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VERIFICATION  OF 
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-  BASE  TEST 
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CURRENT  VELOCITIES 

STATION  i  9 
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TIME,  HOURS  AFTER  MOON'S  TRANSIT  OF  I24TH  MERIDIAN 

MODFL  TEST  DATA 


TIOE  MEAN 

FRESH-WATER  DISCHARGE  11,400  CFS 
SOURCE  fALINIT\  330  PPT 
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-  BASE  TEST 
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CURRENT  VELOCITIES 

STATION  2 1 
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EBB  FLOOD  E~3  FLOOD  EBB  FLOOD 

VELOCITY  FPS  (PROTOTYPE) 
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EBB  FLOOD  EBB  FLOOD  EBB  FLOOD 
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FRESH-WATER  DISCHARGE  11,400  CFS 
SOURCE  SALINITY  330  PPT 


BASE  TEST 


BASE  TEST 
CURRENT  VELOCITIES 
STATION  26 
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EBB  FLOOD  EBB  Fi-OOD  ebb  flood 
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VELOCITY,  FPS  (PROTOTYPE) 


TIDE  .  .  . 

FRESH-WATER  DISCHARGE 
SOURCE  SALINITY. 
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MODEL  TEST  DATA 

TI°E  .  . MEAN  BASE  TEST 
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SOURCE  SALINITY  ...  33.0  PPT  CURRENT  VELOCITIES 
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TIOt  .  MEAN 
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SOURCE  SALINITY  33.0  PPT 

LEGEND 

-  BASE  TEST 


BASE  TEST 
CURRENT  VELOCITIES 

STATION  36 
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VELOCITY,  FPS  (PROTOTYPE) 
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MOPCL  TEST  DATA 


TIDE  MEAN 
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SOURCE  SALINITY  330  PPT 
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BASE  TEST 
CURRENT  VELOCITIES 
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MODEL  TEST  DATA 

TIDE  . MEAN 

FRESH-WATER  DISCHARGE  ...  11,400  CFS 
SOURCE  SALINITY .  330  PPT 
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CURRENT  VELOCITIES 
STATION  4  4 
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PLATE  142 


EBB  FLOOD  EBB  FLOOD  EBB  FLOOD 

VELOCITY,  FPS  (PROTOTYPE) 


VELOCITY,  FPS  (PROTOTYPE) 

EBB  FLOOO  EBB  FLOOD  EBB  FLOOD 


2o7 


PLATE  143 


ebb  FLOOO  EBB  FLOOD  EBB  FLOOD 

VELOCITY,  FPS  (PROTOTYPE) 


VELOCITY,  fps  (prototype) 

ebb  FLOOD  EBB  Fi_OOD  ebb  flood 
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tide  mean 
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SOURCE  SALINITY  33.0  PPT 

LEGEND 
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CURRENT  VELOCITIES 
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SOURCE  SALINITY  33  0  PPT 
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TIME,  HCURS  AFTER  MOON'S  TRANSIT  OF  124°  MERIDIAN 


MODEL  TEST  DATA 

TIDE  MEAN 

FRESH-WATER  DISCHARC-E  11,400  CfS 
SOURCE  SALINITY  33.0  PPT 

LEGEND 

-  BASE  TEST 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATION  40 

SURFACE,  MIDDtPTH  AND  BOTTOM 
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PLATE  180 


TOTAL  SALT,  PARTS  PKR  THOUSAND 


MODEL  TEST  DATA 

T!Ot  .  .  .  MEAN 

FRESH-WAfER  OISCHARGC  .11,400  CFS 
SOORCt  SALINITY  330  RPT 

LEGEND 

-  6ASE  TEST 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATION  44 

SURFACE,  MIDDEPTH  AND  BOTTOM 


PLATE  183 


TOTAL  SALT,  PARTS  PER  THOUSAND 


SURFACE 


TIME,  HOURS  AFTER  MOON'S  TRANSIT  OF  124°  MERIDIAN 


MODEL  TEST  OATA 

TIDE  .  MEAN 

FRESH-WATER  OISCHARCF  ...11,400  CFS 
SOURCE  SAulNITY..  .  .  33.0  PPT 

LEGEND 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATION  45 

SURFACE,  MIOOEPTH  AND  BOTTOM 


BASE  TEST 


PLATE  185 


TOTAL  SALT.  PARTS  PER  THOUSAND 


SURFACE 


MIDDEPTH 


TIME,  HOURS  AFTER  MOON'S  TRANSIT  OF  124 6  MERIDIAN 


MQPEL  TE$T  DATA 
TIOE  MEAN 

fresh-water  oischarge  ii.aoo  cfs 

SOURCE  SALINITY  33  0  PPT 

LEGEND 

-  BASE  TEST 
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SALINITY  OBSERVATIONS 

STATION  47 

SURFACE,  MIDDEPTH  AND  BOTTOM 
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TIME,  HOURS  AFTER  MOON’S  TRANSIT  OF  124°  MERIDIAN 


MODEL  TEST  DATA 

TlOE  MEAN 

FRESH-WATER  DISCHARGE  .11,400  CFS 


SOURCE  SALINITY 


330  PPT 


LEGEND 
-  BASE  TEST 
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SALINITY  OBSERVATIONS 

STATION  48 

SURFACE,  MIDDEPTH  AND  BOTTOM 
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PLATE  187 


TOTAL  SALT,  PARTS  PER  THOUSAND 
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MODEL  1  EST  DATA 
r:OE  .  MEAN 

fhesh-watcr  piSchahce  ;i,aoocfs 

SOURCE  SALINITY  .  33.C  PPT 

LEGEND 

-  BASE  TEST 
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SALINITY  OBSERVATIONS 

STATION  49  . 

SURFACE,  MIDDEPTH  AND  BOTTOM. 


PLATE  188 


TOTAL  SALT,  PARTS  PER  THOUSANO 


40 


SURFACE 


TIME,  HOURS  AFTER  MOON'S  TRANSIT  OF  I.340  MERIDIAN 


MODEL  TEST  DATA 
TIDE  MEAN 

'fresh-water  discharge  H.400  crs 

SOURCE  SALINITY  330  PPT 

LEGEND 

-  BASE  TEST 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATION  51 

SURFACE,  MIDDEPTH  AND  BOTTOM 


PLATE  190 


BOTTOM 


MC-^EL  TEST  DATA 
TIDE  MEAN 

FRESH-WATER  DISCHARGE.  .11,400  CFS 
SOURCE  SALINITY  330  PPT 

LEGEND 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATION  52 

SURFACE,  MIDDEPTH  AND  BOTTOM 


BASE  TEST 


TOTAL  salt,  parts  per  thousand 


40 


SURFACE 


MIDDEPTH 


BOTTOM 


TIME,  HOURS  AFTER  MOON'S  TRANSIT  OF  (£4 0  MERIDIAN 


MODEL  TEST  DATA 

TIDE  MEAN 

FRESH-WATER  DISCHARGE  11,400  CFS 
SOURCE  SALINITY  33  0  PPT 

LEOEND 

-  BASE  tfst 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATION  54 

SURFACE,  MIDDEPTH  AND  BOTTOM 


TOTAL  SALT,  PARTS  PER  THOUSAND 


SURFACE 


BOTTOM 


MODEL  TEST  DATA 

TIDE  MEAN 

FRESH-WATER  DISCHARGE  Il.dOOCFS 
SOURCE  SALINITY  33  0  PPT 

LEGEND 

-  BASE  TEST 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATIC  N  55 

SURFACE,  MIODEPTH  AND  BOTTOM 
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TOTAL  SALT,  PARTS  PER  THOUSANO 
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TIME  *  HOURS  AFTER  MOON'S  TRANSIT  OF  !24°  MERIDIAN 


MODEL  TEST  DATA 

TIDE  MEAN 

FRESH-WATER  DISCHARGE  .11,400  CFS 
SOURCE  SALINITY  330  PPT 


BASE  TEST 


BASE  TEST 

SALINITY  OBSERVATIONS 

STATION  58 

SURFACE,  MiDDEPTH  AND  BOTTOM 
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BASE  TEST  STATION  56 
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